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Abstract 

Within the framework of a two-band tight-binding model, we have performed calcula- 
tions of giant magnetoresistance, exchange coupling and thermoelectric power (TEP) for 
a system consisting of three magnetic layers separated by two non-magnetic spacers with 
the first two magnetic layers strongly antiferromagnetically exchange-coupled. We have 
shown how does the GMR relate with the corresponding regions of magnetic structure 
phase diagrams and computed some relevant hysteresis loops, too. The GMR may take 
negative values for specific layers thicknesses, and the TEP reveals quite pronounced 
oscillations around a negative bias. 
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1 Introduction 



Since the discovery of giant magnetoresistance (GMR) phenomenon in magnetic multilayers |l] 
there has been a great deal of interest in studying them by both theoretical and experimental 
methods. The reason for the interest is, already partially realized, possibility of practical 
applications as magnetic sensors, recording heads and magnetic memory elements. 

The standard system exhibiting the GMR is a trilayer (i.e. two magnetic layers separated 
by a non-magnetic spacer) with thickness of the spacer chosen so as to produce antiferro- 
magnetic coupling between magnetic layers. While such a system, due to its simplicity, is 
convenient for theoretical treatment, it presents some problems in practical applications. The 
main problem is a high switching field which is usually necessary to rotate the magnetizations 
(overcoming antiferromagnetic coupling) and produce GMR. One way to deal with this prob- 
lem is to use somewhat more complex spin- engineered structures. Widely known structures 
of this type are spin- valve systems 0, in which one of the magnetic moments is fixed by the 
strong exchange coupling due to an additional antiferromagnetic layer (e.g. MnFe or CoO). 
There exists however also a different approach in which a system composed of three magnetic 
layers is used. Two of them are strongly antiferromagnetically coupled, forming the so called 
artificial antiferromagnetic subsystem (AAF) ||, and the third one — detection layer — is 
only weakly coupled (or just decoupled). Such a setup was proposed both for laboratory 
measurements @, || and more practically as angular velocity meter ||. A similar system 
(superlattice with strong and weak exchange couplings) was also studied theoretically (on ab 
initio level) 0, the thicknesses involved were however small due to numerical limitations. 

The aim of the present paper is to perform thorough studies of transport and magnetic 
properties of the systems in question and to relate them with corresponding magnetic structure 
phase diagrams. 

2 The model and the method of calculations 

We consider a system consisting of three magnetic layers separated by two non-magnetic 
spacers, i.e. the structure of the F\/ 'Si/ 'F 2 / Sz/ 'F 3 type, where Fi stands for ferro- and Si 
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for paramagnetic layer. In order to describe collinear configurations we employ tight-binding 
hamiltonian with two, hybridized bands and spin-dependent on-site potentials (see Ref. JF| 
for details). We restrict ourselves for simplicity to the case of simple cubic structure. The 
following values of model parameters have been chosen: Ep = 0, t s = —1, = —0.2, V s d = 1, 
t\ a = 0, ef| = —1 for all the layers, ef^ = —0.2 within magnetic layers and —1 elsewhere, 
where Ep is the Fermi energy, t a (with a being the band index — s, d) are the hopping 
integrals, V s d is the s — d intra-atomic hybridization and e" CT are the on-site potentials (where 
a is \ for majority- and J. for minority-spin carriers). The above set of parameters enables 
us to mimic essential features of the electronic structure of Co/Cu multilayers {i.e. the spin- 
polarized density of states is qualitatively reproduced — in particular the majority <i-bands 
in magnetic and (non-polarized) rf-bands in paramagnetic layers are matched perfectly which 
closely resembles the situation in Co/Cu systems). 

The conductance is computed from the Kubo formula with the help of recursion Green 
function technique |7|, |8|. The only difference in comparison with |7j is that hybridization in 
lead wires (attached to the multilayer from both sides, for transport calculations) has been 
taken into account this time. The GMR has been defined as 

rU;J. 

GMR =nF- 1 ' 

where the arrows show the orientations of magnetic moments. Note that without the first mag- 
netic layer this definition would be identical with the usual one. Additionally the thermopower 
(TEP) has also been calculated from the following formula (see eg. Ref. ||) 

We define, as in Ref. [0, the "giant magneto-TEP-effect" GMTEP analogously to Eq. (0) 
(with T replaced by S). 

For studying the magnetization processes we employ the phenomenological expression, not 
unlike the one introduced in Ref. ]TTJ, (Gj is an angle between the z-th magnetic moment and 
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the external field, B) 



E (61,62,63) = -J12 COS (61-62)- J 23 COS (62-63)- Jl3 008(0! -63) 

3 3 

-B J2 U cos (0 4 ) /t + J2 £a(6;), (3) 

1=1 4 = 1 

where the first three terms describe bilinear exchange coupling between magnetic layers, the 
next three are Zeeman energy terms (ij being z-th layer thickness and t the overall thickness of 
all the magnetic layers) and Ea{Qi) is the crystalline anisotropy, that is tiK sin 2 (26j)/4t and 
— tiD cos 2 Qi/t for cubic and uniaxial case respectively. We assume that external magnetic 
field is applied along the [10] in-plane crystallographic axis, which can be either easy or hard 
axis depending on the sign of the anisotropy constants. The magnetic moments are confined 
to the layers plane which corresponds to the strong shape anisotropy. Expression (|3|) was then 
numerically minimized, with respect to the 6j-s, by taking, starting from initial configuration, 
little steps in the direction opposite to the energy gradient. All the extremal points found in 
this way were additionally checked against the stability condition (i.e. the positivity of all the 
minors of Mjj = d 2 E/(dOidOj)) in order to eliminate saddle points. Using Eq. @ one can 
write (for B equal to 0) 

J12 = - [E (ttOO) + E (OttO) - E (000) - E (OOtt)] , 

J 23 = - [E(00n) + E (OttO) - E (000) - E (vrOO)] , (4) 

Ji 3 = \\ E ( 007r ) + E (too) - E (00°) - E ( 07r0 )] • 

Therefore, having known the energies of collinear configurations from the model calculations 
(based on the two-band tight-binding hamiltonian) , we are able to determine the exchange 
coupling constants. 

From now on we will be using reduced values of the magnetic field b = -B / [ «7i2 1 and the 
anisotropy constants k = K/\ J 12 | and d = D/\ «/ 12 l ■ We will also assume, if not stated explicitly 
otherwise, the magnetic layers thicknesses to be 8, 3 and 3 ML (monolayers), respectively, in 
order to keep the length ratios as in Ref. ||. The first spacer thickness will be set to 3 ML in 
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order to achieve the needed antiferromagnetic coupling between the first two magnetic layers. 



3 Results 

Figure [3] presents the exchange coupling constants ( J-s) plotted against the thickness of the 
second spacer (ns2)- As mentioned above, for the chosen thicknesses we got strong antifer- 
romagnetic coupling (J12) between the first two magnetic layers, while J23 and J13 oscillate 
around zero. In all three cases the period of oscillations is about 3 ML, which is close to the 



theoretically predicted value (2.8 ML) coming from the stationary spanning vector |I3 
placed at the (0, ir) (and equivalent positions) in the two-dimensional Brillouin zone. The sec- 
ond period (8 ML) originating from the hole pocket placed at (71", ir) does not seem to appear 
in the present context. Note however that it can become visible under some circumstances 
like in Ref. [0 where the tunnelling conductance has been considered. The GMR and J23 for 
the same system have been plotted in Fig. 0. It can be noted that GMR asymptotically tends 
to oscillate with the same period but in opposite phase to J23- This is in agreement with 



our previous findings fl3| , but it is still not clear to what extent this correlation is universal. 
The values of GMR are strongly reduced in comparison with the trilayer case. This can be 
easily understood if we note that, due to the fixed antiferromagnetic alignment of the first 
two magnetic layers, there is no non-scattering channel in any of the configurations involved 
(see Eq. [1|). Figure |3] where the on-site potentials for the c?-bands (ef a ) have been schemat- 
ically plotted, shows that there exist at least two scattering interfaces in each case. Basing 
on the number of interfaces one can qualitatively predict that the f |; J, down- and ||; | up- 
spin electron channels have the higher conductances than the remaining two. This is indeed 
clearly visible in Fig. |] where the computed conductances are shown. As already discussed 
there is no obvious highest conductance channel. Instead, we have two higher and two lower 
conducting channels close to each other within the pairs. As a consequence the sign of GMR 
is determined by all the channels, and can be changed by manipulating some parameters (eg. 
thicknesses of the layers). This is the case in Fig. |5] where we have plotted the GMR and 
J23 for a system with thickness of the second magnetic layer set to 5 ML (instead of 3 ML 
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as in Fig. [|). The GMR oscillations have a small but clearly negative bias. The asymptotic 
opposite-phase correlation with respect to J 2 3 is again clear in this case. 

The GMTEP, calculated for the same set of parameters as in Fig. [| has been plotted in 
Fig. []. In agreement with the findings of Ref . [[Uj] , the oscillations are quite pronounced and 
have the same period as GMR but exhibit negative bias. Asymptotically they seem to have 
roughly the same phase as GMR. 

For studying the magnetization processes we have chosen the parameters as in Fig. ^ with 
the second spacer thickness (ns^) set to 7 ML (the second ferromagnetic maximum of ^23)- As 
already stated, two cases have been taken into account, i.e. the cubic and uniaxial anisotropies. 

In the first case the anisotropy term (tj/csin 2 (2Bj)/4t) gives rise to four potential wells 
placed at the following in-plane crystallographic axes : [10], [01], [10] and [01] for k > and 
[11], [ll], [II] and [11] for k < 0. Figure [7|a exemplifies some phase diagrams for various initial 
configurations. We have chosen the simplest ones i. e. the collinear configurations with relative 
alignment of magnetic moments favoured by interlayer exchange coupling. It is however, in 
principal, possible to stabilize also non-collinear ones, provided that the anisotropy is strong 
enough. The phase diagrams have been obtained by taking subsequent scans along b for 
different k values. Dotted horizontal lines are thus only guides to the eye. The diagrams exhibit 
rich structure with a number of configurations (phases) occurring during the magnetization 
process (including the non-collinear configurations). The transition between them can be either 
of first or second type, that is they manifest themselves as discontinuities in magnetization 
or its first derivative. Some exemplary hysteresis loops are presented in Fig. [7|b. As expected 
stronger anisotropy produces a richer structure. For positive, and sufficiently big, values of k 
some flat regions, typical for exchange-biased systems, occur. Note however that there is no 
11] I t° Hi T transition in the first diagram of Fig. 0a (but see below). 

For the case of uniaxial anisotropy (— t id cos 2 (<di)/t) there exist only two potential wells, 
that is [10] and [10] for d > and [01], [01] for d < 0. The phase diagrams (Fig. ||a) are 
somewhat less complicated now, due to the simpler energy landscape. The saturation curve 
for the upper part of the first diagram (with the tli I initial configuration) has been obtained 
on analytical basis (i.e. from the stability condition) because of the weak stability of the |J.; j 
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configuration (by which we mean that the existing energy minimum is very shallow) which 
makes it difficult to perform reliable numerical minimalization. The same precautions have 
been applied to the first hysteresis loop in Fig. [|b. Note that this time the flat regions are 
present already for small values of d. For d > 0.51 there exists, in the first diagram, the above 
mentioned ||; J, to fl; T transition. Only positive values of d were presented since the opposite 
case is trivial — there is practically no hysteresis. 

4 Conclusions 

Within the microscopic two-band tight-binding model we have performed the calculations of 
interlayer exchange coupling, current-perpendicular-to-plane conductance and thermopower 
for a system consisting of three magnetic layers, separated by paramagnetic ones. With 
the thicknesses chosen so as to produce strong antiferromagnetic coupling between the first 
two magnetic layers, we found both the interlayer exchange coupling, GMR and GMTEP to 
oscillate, as a function of the second spacer thickness, with the period originating from one of 
the extremal spanning vectors of the Fermi surface. Additionally, using the phenomenological 
approach, we have computed magnetic phase diagrams and commented on their relevance to 
the magnetoresistance. We found that the phase diagrams exhibit rich structure and there are 
flat regions in hysteresis loops, typical for exchange-biased spin-valves. It has been also found 
that in the case of the systems under consideration, in contrast to conventional trilayers, it is 
possible to obtain a negative (inverse) perpendicular GMR by merely changing thicknesses of 
particular layers. 
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List of Figures 

"Exchange coupling constants, calculated for the Tj / Si / F 2 / S "2 / F 3 structure" 



Exchange coupling constants, calculated tor the b 1/01/^2/02/ -^3 structure 
plotted against the second spacer thickness (ns 2 )- Magnetic layer thicknesse ; 



have been set to 8, 3 and 3 ML, respectively and the first spacer thickness is 



equal to 3 ML. The hamiltonian parameters are as follows: Ep = (Fermi en- 



ergy), t s = — 1, td = —0.2 (hopping integrals), V s d = 1 (hybridization), e s in ="0 



and eft = — l(for all the layers), efj = —0.2 for i within the magnetic layer and 



-1 elsewhere. Note the strong antiferromagnetic coupling (Ji 2 ) between the 



first two magnetic layers.] 11 



2 CPP-GMR and the interlayer exchange coupling (J23) plotted as a function 



of the second spacer thickness. The period of oscillation is about 3 ML in 



agreement with the predicted value of 2.8 ML, resulting from the corresponding 



Fermi surface calliper. Asymptotically the oscillations of both quantities are 



predominantly opposite in phase.] 12 



3 Schematic plot of the on-site potentials for d-bands (efj across the b\j ' S\j ' b\j ' S2I b% 



multilayer for different configurations and spin directions. The tli I down- and 



\l; t up-spin channels with only two scattering interfaces are clear candidates 



for high conductance channels. Note that there is however no non-scattering 



channel for configurations taken into account .| 13 



4 Conductances (in e A /h units) for parameters as in Fig. 2\ There exist, as dis- 



cussed in the text, two high and two low conductance channels 13 

5 CPP-GMR and the interlayer exchange coupling ( J23) calculated for the same 



set of parameters as in Fig. 2| , except for the thickness of the second magnetic 



layer (n/2) which has been set to 5 ML. The GMR oscillate around a small 

negative value 14 

5 GMTEP (giant magneto-TEP effect) for the same set of parameters as in Fig. % 



GMR has been also plotted for reference. The GMTEP oscillates asymptotically 

in roughly the same phase but around a negative bias 15 
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7 a) Phase diagrams for magnetization processes obtained for different initial con^ 
figurations in the case of cubic anisotropy. The second spacer thickness (rasJJ 
has been set to 7 ML. The arrangements of magnetic moments are schematic 



cally drawn with the first (8 ML) magnetic layer symbolized by the long solid 
arrow and the third one (detection layer) by the dotted arrow, b) Exemplary 



hysteresis loops with the dashed line marking the backward sweep. The zero 



field configurations are schematically marked for both the sweep directions. 



For positive values of fc, the strong anisotropy produces flat regions typical for 



exchange-biased systems. 

8 a) Phase diagrams in the uniaxial anisotropy case, b) Exemplary hysteresis 



loops. Only positive values of d are shown since the opposite case exhibits no 
hysteresis at all. The flat regions are present already for weak anisotropy. . . . 
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Figure 1: Exchange coupling constants, calculated for the F x j S\j F 2 / S 2 / F 3 structure, plotted 
against the second spacer thickness {ns2}- Magnetic layer thicknesses have been set to 8, 
3 and 3 ML, respectively and the first spacer thickness is equal to 3 ML. The hamiltonian 
parameters are as follows: Ep = (Fermi energy), t s = —1, td = —0.2 (hopping integrals), 
V sd = 1 (hybridization), e| CT = and ef T = — l(for all the layers), efy = —0.2 for i within the 
magnetic layer and —1 elsewhere. Note the strong antiferromagnetic coupling (J12) between 
the first two magnetic layers. 
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Figure 2: CPP-GMR and the interlayer exchange coupling (J23) plotted as a function of the 
second spacer thickness. The period of oscillation is about 3 ML in agreement with the pre- 
dicted value of 2.8 ML, resulting from the corresponding Fermi surface calliper. Asymptotically 
the oscillations of both quantities are predominantly opposite in phase. 
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Figure 3: Schematic plot of the on-site potentials for d-bands (ef a ) across the F 1 / Si/ F 2 / S 2 / F 3 
multilayer for different configurations and spin directions. The ||; J, down- and ||; | up-spin 
channels with only two scattering interfaces are clear candidates for high conductance channels. 
Note that there is however no non-scattering channel for configurations taken into account. 
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Figure 4: Conductances (in e 2 //i units) for parameters as in Fig. ||]. There exist, as discussed 
in the text, two high and two low conductance channels. 
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Figure 5: CPP-GMR and the interlayer exchange coupling ( J23) calculated for the same set 
of parameters as in Fig. |2], except for the thickness of the second magnetic layer (n/2) which 
has been set to 5 ML. The GMR oscillate around a small negative value. 



14 



0.024 




0.016 



Figure 6: GMTEP (giant magneto-TEP effect) for the same set of parameters as in Fig. |[ 
GMR has been also plotted for reference. The GMTEP oscillates asymptotically in roughly 
the same phase but around a negative bias. 
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Figure 7: a) Phase diagrams for magnetization processes obtained for different initial config- 
urations in the case of cubic anisotropy. The second spacer thickness (ns2) has been set to 
7 ML. The arrangements of magnetic moments are schematically drawn with the first (8 ML) 
magnetic layer symbolized by the long solid arrow and the third one (detection layer) by 
the dotted arrow, b) Exemplary hysteresis loops with the dashed line marking the backward 
sweep. The zero field configurations are schematically marked for both the sweep directions. 
For positive values of k, the strong anisotropy produces flat regions typical for exchange-biased 
systems. 
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Figure 8: a) Phase diagrams in the uniaxial anisotropy case, b) Exemplary hysteresis loops. 
Only positive values of d are shown since the opposite case exhibits no hysteresis at all. The 
flat regions are present already for weak anisotropy. 
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